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proper t loo . The rosonrch of the Saharan dust will lend to a 
siBiilf leant study in the future when the mode 1 will be avollnblo to 
couple with nlr chemistry, nnd cloud and precipitation physics. This 
complete aerosol model would explain hopefully the micro- and macro- 
physical properties result liter from the complete cycle of aerosol 
evolutions, nnd the nna lysis of the validated fields against 
observations would help to estimate the effects on clouds, on 
radiative clutrae ter is t len (Patterson, Cillottc and Stockton, 1977: 
Carlson nnd Cover ly, 1977), and on visibility. 
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THEORETICAL STUDIES OF THE TRANSPORT AND REMOVAL PROCESSES 

OF SAHARA DUST 


A subsynopt lc scale aerosol model hao boon dve loped during- the 
Inst six month at Ames Research Center to study the long term fate of 
the Snliarnn dust which may piny an important role In affecting 
the climate over the Atlantic and Europe by cooling the atmosphere 
due to the increase of the earth’s albedo or conversely warming the 
atmosphere due to the enhancement of the greenhouse effect. 

I 

The processes controlling the evolution of aerosol particles may 
bo divided into dynamic processes concerned with the motions of air 
currents and microphyo lea 1 processes concerned with the growth of the 
individual particles. If dust particles are introduced Into the 
atmosphere, the particles will be diffused through a certain layer 
and transported along the air currents. At the 800 mb level over the 
Atlantic (sec Figure 1. Gates, 1970) the dust particles originating 
from North Africa will be transported west at inland, southwest near 
African coast and then be transported northwest passing the ridge. 
However, the situation is quite different at the 400 mb level (see 
Figure 2. Gates, 1976) at which the particles whose concentration may 
be quite small will bo transported along the westerly jet to Europe 
and Mid-East. 

The main body of the dust clouds may be contained In the 
convective layer and transported east to reach the eastern coast of 
North and Central America. The rate of production of the large 
particles depends on the collision efficiency which depends in turn 
cn the flow of air In which the particles are Imbedded and the 
dynamics cf the particles in response to the drag force exerted by 



the air (Langmuir, 1948; Penrcey and Ilill, 1957; Hocking, 1959 ; Shafrir 
nnd Ho lburger , 1963 ; Davis and Sartor, 1967; Davis, 1972) as well as 
tlielr coalescence efficiency. Studies have been made in order to 
parouie tor Ize the removal processes due to coagulation and 
sod imontot ion, 

A transport jtiodo 1 lias been developed based on th" Ames Planetary 
Boundary Layer Model (Loo and Swan, 1977; Swan and Lee, 1978). The PBL 
model has been modified to Include dust transport in the convective 
layer. Therefore, the model consists of three layers; mixed layer, 
convective layor nnd free atmosphere. The dust primarily blows off 
the African desort in a vresterly direction and the dust lnyer 
thickness secm3 to be controlled by the daytime depth pf the 
planetary boundary layer over Africa, and remains constant over the 
ocenn. By the addition of n third layer to the model, this 
convective layer transport can be taken into account. 

The results from a preliminary tost run arc presented In 
Figures 3 to 6. For this test run, simple initial and boundary 
conditions were assumed to provide PBL wind n3 southeasterly with 6 
m/scc, sigma value at interface os 0.85 (about 1200 meters above 
surfnee), PBL and surface potential temperature as 295 K, potential 
temperature Jump at interface as 10 K, stability In the convective 
layer as 50 degrees per sigma and vertical gradients o f x- and y- 
component wind as -40 and 40 m/sec- per sigma, respectively. Those 
simulated fields which represent the subsynoptic scale motions for 
4!00 a.m. of July 15 were obtained after 48 hours of real tlmb 
simulation In order to pass the transient period. The surface 
pressure field (Figure 3) shows the topographic features of the 
target area as we 1 1 as the mass distribution of air due to the 
initial southeasterly wind. The vertically averaged parameters of PBL 
depth, potential temperature and wind are shown in Figures 4 to 6. 

The wind pattern over ocean was quite uniform with about 5 m/sec, but 


the wind over land was channelled, blccTod aud drained due mainly lo 
tlie topographic features. The actual simulation of those paramo tops 
t/ 111 be quite different If the Initial nnd boundary valuos that 
varies with time nnd space are introduced. 

In order to simulate tlie real situation as closely na possible, 
the Initial values cud tlie steady state boundary values are provided 
from various sources of observational and numerical studies (Posey 
and Clapp, 1904; Espensliadc, Jr., 1970; Alexander and Hobley, 1970; 
Gates, 1976) and these row data wore interpolated and smoothed to ftt 
properly to tlie coordinate system employed by the model nnd to tlie 
complex terrain of the target area, in Figures 7 nnd 0, the surface 
albedo constructed by Posey and Clapp( 1964) and the surface heat 
capacity field constructed accounting the soli properties 
( Espensliadc , 1970) over the target area are shown, respectively. In 

Figures 9 through 22, tlie self-explanatory Initial fields constructed 
to fit tlie model ore shown (Gates, 197G; Alexander and Kobley, 1976) 
for various parameters. The actual simulations using these initial 
values are not available at tills time. However, we can discuss a few 
interesting features found from them. As expected, the hot and dry 
spo tG (Figures 10 ami 10) are located over most of Africa. Tlie PBL is 
shallow over loud. In particular It Is shallower over Saharan 
desert. It Is quite moist In the PBL (Figure 13) and qulto dry In 
the convective layer (Figure 14). .The PEL momentum field (Figures 15 
aud 16) can be explained considering the surface pressure pattern 
(Figure 9) . 

In the future, this model may be exercised extensively to 
understand various dynamic and physical processes characterising the 
behavior of the dust and to validate the results against 
observations. The surface dust Inventory, If available, will be used 
to parameterize tlie dust flux In terms of mesoscale parameters such 
ns momentum, drag and atmospheric stability and of surface 
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Figure 5. PBL Potential Temperature in K (Test Run) 
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Figure 6. PBL Wind with Maximum 10 m/sec (Test Run) 
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Figure 11. Sigma Value at 
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